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\Q • ABSTRACT 

in , 

Aims. The progenitor stars of long Gamma-Ray Bursts (GRBs) are thought to be Wolf-Rayet stars, which generate a massive and energetic 
so . wind. Nevertheless, about 25 percent of all GRB afterglows light curves indicate a constant density medium close to the exploding star. We 

explore various ways to produce this, by creating situations where the wind termination shock arrives very close to the star, as the shocked wind 
^ C"| material has a nearly constant density. 

O | Methods. Typically, the distance between a Wolf-Rayet star and the wind termination shock is too large to allow afterglow formation in the 
shocked wind material. Here, we investigate possible causes allowing for a smaller distance: A high density or a high pressure in the surrounding 
interstellar medium (ISM), a weak Wolf-Rayet star wind, the presence of a binary companion, and fast motion of the Wolf-Rayet star relative 
to the ISM. 

, Results. We find that all four scenarios are possible in a limited parameter space, but that none of them is by itself likely to explain the large 
• • ■ fraction of constant density afterglows. 
, Conclusions. A low GRB progenitor metallicity, and a high GRB energy make the occurrence of a GRB afterglow in a constant density medium 
, more likely. This may be consistent with constant densities being preferentially found for energetic, high redshift GRBs. 

C« Key words. - Stars: Wolf-Rayet - Stars: winds, outflows - Gamma rays: bursts - ISM: bubbles - Hydrodynamics - 
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1. Introduction 

During their evolution, massive stars lose a major fraction of 
their mass in the form of a stellar wind. The interaction between 
this stellar wind and the surrounding interstellar medium cre- 
ates a circumstellar bubble. Within approximations, the mor- 
phology and evolution of such bubbles can be predicted ana- 
lytically (Castor et al. I1975I Weaver et al. I1977I Ostriker & 
McKee 1988 1, or by numerical simulations (Garcia-Segura et 
al. I1996al I1996bl Freyer et al. E0TB1 van Marie et al. 2004 
2005a 2005b I. Figure^shows schematically the morphology 
of a spherical wind bubble. Closest to the star is the free- 
streaming stellar wind. Its material hits the wind termination 
shock, where the kinetic energy of the wind material is ther- 
malized, which feeds the 'hot bubble' of shocked wind mate- 
rial. Outside the hot bubble lies a shell of shocked interstellar 
gas that has passed through the outer shock that delimits the 
bubble as long as it expands supersonic ally with respect to the 
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interstellar medium (ISM). The material forming this shell is 
swept up as the high-pressure gas of the hot bubble. Finally, 
there is the ambient medium, which has not yet been overrun 
by the expanding circumstellar bubble. 

Within the collapsar picture (Woosley 119931 . Wolf-Rayet 
stars, which form the final evolutionary stage of massive stars, 
are thought to be the progenitors of long Gamma-Ray Bursts 
(GRBs). This relativistic jet associated with the GRB expands 
into the circumstellar medium (CSM). On its way it first en- 
counters the free-streaming wind of the Wolf-Rayet star, and 
only afterwards it may penetrate into the shocked wind mate- 
rial. 

The medium close to a Wolf-Rayet star has the density pro- 
file of a free-streaming wind with a total mass flux M and ve- 
locity V w : 



The wind density drops as p w oc r 2 once the wind has reached 
its terminal velocity, which occurs at a distance of a few stellar 
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radii. Both analytical and numerical calculations suggest that 
the free wind of a Wolf-Rayet star typically extends over a ra- 
dius of the order of a few parsec. This is confirmed by observa- 
tions of Wolf-Rayet nebulae, such as NGC 6888 and RCW 58, 
which also have radii of the order of several pc (Gruendl et al. 
120001 . 

However, from the analysis of many GRB afterglows, a 
constant circumstellar medium density is inferred. The fraction 
of GRBs that show a constant density profile varies, depending 
on the afterglow model, but it can be estimated to > 25 percent 
(Chevalier & Li 2000 Panaitescu & Kumar 2001 Panaitescu & 
Kumar. 120021 Chevalier et al. 2004 1. Moreover, it is notewor- 
thy that all studies mentioned were limited to GRBs detected 
before the launch of the SWIFT satellite, i.e. they are primar- 
ily concerned with low redshift GRBs. It cannot be excluded 
that the fraction of afterglows in a constant density medium is 
higher for larger redshift. 

The density of the CSM and its radial dependence is also 
relevant for the use of GRBs as standard candles in cosmo- 
logical studies. For this purpose, the correlation between peak 
energy and corrected collimated energy has to be calculated, 
which depends on the nature of the circumstellar or interstellar 
medium (Ghirlanda et al. 2005 1. 

Here, we investigate various possibilities to obtain GRB af- 
terglows in a constant density medium. In Sect. [2]we estimate 
analytically at what radius the constant density medium needs 
to occur in order to appear in the GRB afterglow lightcurve. 
In Sects. [3]and|4] we explore ways to produce a constant den- 
sity medium sufficiently close to the star, using both analytical 
and numerical techniques to model the circumstellar medium. 
In Sect.[5]we discuss the likelihood of the various possibilities. 

2. Analytic considerations 

While the region closest to the Wolf-Rayet star has a density 
that decreases with the radius squared, the next layer of the 
circumstellar bubble (region B in Fig. Q has a nearly con- 
stant density, essentially since the high sound speed inside 
this material precludes the occurrence of large pressure- and 
density-gradients. This hot bubble of shocked wind gas drives 
a shell into the CSM (Weaver et al. U977l l. It has been proposed 
(Wijers et al. l200Tl Chevalier et al. l2004l that the GRB after- 
glows are created while the jet moves through this medium. 
This requires the Gamma-Ray Burst jet to cross the wind ter- 
mination shock at an early stage. Since Wolf-Rayet stars have 
a strong wind, with a high mass loss rate and a high termi- 
nal velocity, the ram pressure of the wind is considerable, and 
the wind termination shock that separates this bubble from the 
free-streaming wind is typically located at a distance of 10 pc 
from the star, which would make this scenario impossible. (For 
a discussion of the models by Chevalier et al. (2004i, which 
find considerably smaller shock radii, see Sect. 3.3). 

2.1. The location of the wind termination shock 

The location of the wind termination shock can be calcu- 
lated analytically, using the so-called thin-shell approximation 
(Castor et al. \Wl5[ Weaver et al. [T977I Ostriker & McKee 
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Fig. 1. Schematic view of a circumstellar bubble (not to scale). 
The free-streaming stellar wind (region A) passes through the 
wind termination shock R w (dashed line) to enter the hot bub- 
ble of shocked wind material (region B). The high thermal pres- 
sure in the hot bubble sweeps up a shell (thick, continuous line) 
of radius R\,, which expands into the ISM (region C). 



1988 1. This assumes that the forward shock is completely ra- 
diative, whereas the hot bubble of shocked wind material is 
adiabatic, which serves as a good first approximation. The me- 
chanical luminosity of the star, 

K = - M w Vl , (2) 

is a measure of the mechanical energy that the star puts into the 
CSM. For a given mechanical luminosity, it is possible to cal- 
culate the thermal pressure in the bubble as a function of time, 
by balancing the pressure in the hot, shocked wind with the 
ram pressure po(d/?b/df) 2 due to the expansion into a constant- 
density ISM with mass density po (Weaver et al. U977l l: 

Pb = Z_ L 2/5 3/5 r 4/5^ (3) 

(3850tt) 2 / 5 w h ° 

This ansatz assumes that the ISM is cold, so that it has no sig- 
nificant thermal pressure. Since the mechanical luminosity of 
the star can change considerably over time, we have to take the 
average value for the mechanical luminosity. 

The wind termination shock is located at the radius R w such 
that the thermal pressure in the shocked wind material equals 
the hot bubble pressure P\,. The post-termination shock pres- 
sure is roughly equal to the ram pressure in the wind, 
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Combining Eqs.[3]and|3]gives the radius of the wind termina- 
tion shock: 



t 2/5 



Mw V w \ 1 / 2 /3850\ 1 / 5 

U' (Wo) 3 ' 10 ' 



(5) 



Here the mass loss rate ( M w ) and wind velocity ( V w ) are time 
dependent, while L w is the average mechanical luminosity over 
the entire evolution of the star up to that moment. To obtain the 
radius of the wind termination at the end of the Wolf-Rayet 
phase, one should take the properties of the Wolf-Rayet wind 
for M w and V w , and one should use the average value of the 
mechanical luminocity. 1 

If the ram pressure of the wind changes, the radius of the 
wind termination shock will change on the hydrodynamical 
time scale of the hot bubble. Since the temperature in the hot 
bubble is of the order of 10 7 K, this timescale is 1...4xl0 4 
years, for typical bubble radii of few tens of parsecs (Garcia- 
Segura et al. U996al[T996bl Freyer et al. !2UU51 van Marie et al. 
2005b). Consequently, the wind termination shock will closely 
follow any changes of the stellar wind parameters which occur 
on the time scale of the stellar evolution. The time dependence 
of the thermal pressure in the bubble is not important here, as 
we only consider the Wolf-Rayet period in the following, which 
lasts only a fraction of the total stellar lifespan. 

2.2. Afterglow production in a free wind 

In order to establish a maximum radius for the wind termina- 
tion shock in those cases where a constant-density circumstel- 
lar medium is inferred, we must first calculate the typical dis- 
tance of the jet from the progenitor star when the afterglow 
is produced. The afterglow of a GRB is believed to be syn- 
chrotron radiation, produced when the jet produces a relativis- 
tic shock in the surrounding medium. For a constant-density 
medium to show up in the afterglow light curve, the jet must 
still be relativistic when it passes through the wind termination 
shock. 

We make an analytical estimate of the distance that the 
GRB jet penetrates into the surrounding medium before slow- 
ing down to non-relativistic speed. C.f. Piran (1999 2005 1, 
Panaitescu & Kumar (2000 1 and Chevalier and Lee (2000 1. 
This provides a rough estimate of the radius at which the con- 
stant density medium must be found. For a more accurate result 
it is necessary to model the GRB afterglow lightcurve gener- 
ated in a structured circumstellar medium. Such models were 
presented by Ramirez-Ruiz et al. (2005 1 and Eldridge et al. 
(2005). 

We consider an extremely powerful fireball with an ex- 
plosion energy £grb — 10 52 ...10 54 erg, and with a small 



1 N.B. This is only true for the situation where the Wolf-Rayet 
wind expands directly into the main sequence bubble. During the early 
Wolf-Rayet stage, the Wolf-Rayet wind sweeps up the wind of the 
preceding phase (Red Supergiant or Luminous Blue Variable). In that 
case, the radius of the wind termination shock can be found by using 
Eqn. (31) from Garcia-Segura and Mac-Low ( 1995 l. 



mass loading in the form of baryons: McRB(baryonic) < 
10 4 ...10 5 M . Note that £ G rb is the energy of the fireball 
for a spherical explosion. The actual energy of the jet can be 
much lower, depending on the opening angle. A jet confined to 
a solid angle f2j carries an energy £j ~ £grb x (Qj/47r) and 
amassMj ~ M G rb X (Oj/4?r). 

As the total explosion energy is significantly larger than the 
rest energy of the ejected mass, 



Egrb » Mgrbc 



(6) 



the material in the blast wave (or jet) will expand rapidly, mov- 
ing with a speed V s « c and with a corresponding bulk Lorentz- 
factor: r s » 1. During the initial phase of the expansion, the 
explosion energy £grb is converted by pressure forces into the 
bulk motion of the material. Once the initial energy released 
by the GRB has been converted into bulk kinetic energy, the 
Lorentz-factor is determined solely by the ratio between the 
explosion energy and the rest-mass energy (see the review by 
Piran, 2005 and references therein): 



n = 



-GRB 



Mgrbc 2 



(7) 



This relation will hold until the jet has swept up a sufficient 
amount of mass (see below), or until (in the case of a jet) side- 
ways expansion becomes important. In order to explain the ob- 
servations, we need 77 =a 100. ..1000 (Thompson 1994 1. 

During the expansion, the blast wave sweeps up material 
which passes through a relativistic shock. The total energy of 
the system, consisting of the original mass and radiation liber- 
ated at the GRB and the swept-up mass, approximately equals 



£grb ~ r s M CRB c 2 + ( r 2 - 1) M sw c 2 



(8) 



The second term is the total (internal) energy of the shocked, 
swept-up mass. For T s » 1 this term follows from the jump 
conditions for an ultra-relativistic shock (Blandford & McKee, 
119761 which give the post-shock energy density as <?2 ~ F 2 pic 2 
where the subscript 1 (2) refers to the pre-shock (post-shock) 
state, and where p\ is the rest mass density in front of the 
shock. For a non-relativistic shock with velocity V s the clas- 
sical Rankine-Hugoniot relations give e2 ~ P1V 2 . The second 
term in expression (|8j correctly interpolates between these two 
limiting cases, up to factors of order unity. £grb is conserved 
as long as radiation losses can be neglected. 
In the limit F s » 1 one can write: 



£grb = T s Mgrbc 2 + T 2 M sw c 2 . 



(9) 



The free expansion phase, where T s ~ £grb / -^grb ° 2 ~ con- 
stant, ends when the energy contained in the shocked, swept-up 
mass becomes comparable with the energy of the original ex- 
plosion: 



Mi 



GRB 



M, 



GRB 



^GRB 
T] 2 C 2 



= M d , 



(10) 



which defines the deceleration mass M&. The blast wave or jet 
will decelerate when M sw > Taking the simple case of a 
spherical blast wave expanding into a freely expanding stellar 
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wind, with a density profile as in Eqn.^ an d constant V w , the 
swept-up mass when the fireball has expanded to a radius R s 
equals 



M W R S 

Vv 



(ID 



Combining relations dlOi and (II H one finds the deceleration 
radius, the radius where the free expansion of the blast wave 
stops, provided this occurs while still in the free-wind region: 



R A = 



M d V w E GRB y w 



(12) 



The blast wave now enters the relativistic equivalent of the 
well-known decelerating Sedov-Taylor expansion, until the ex- 
pansion speed becomes sub-relativistic. After that, it enters the 
ordinary Sedov-Taylor expansion phase until radiation losses 
become important, invalidating the assumption of the conser- 
vation of the total energy (Eqn|HJ. During the Sedov-Taylor 
phase, the swept-up mass dominates the energy equation, so 
that: 



£grb = F 2 M SW c 2 



(13) 



As long as the blast wave stays in the free wind, relation il lb 
remains valid, and one has: 




-1/2 



(14) 



where we have used Eqn.^] We can take the point where F s a 
1 as the point where the blast wave becomes non-relativistic, 
This happens at the so-called Sedov-Taylor radius Rsj, which 
in a free wind equals: 



_ 2 p _ e grb Vy 
«st = rj Rd - — -r — — 

M w c 2 



For R s > R SJ one has from with t 2 - 1 « V 2 /c 2 - 



V 
c 



Rsj 



-1/2 



(15) 



(16) 



for as long as the expansion proceeds in the free-wind region 
(R S <R W ). 

For a typical Wolf-Rayet wind one has: £grb = 5.0 x 
10 52 erg, V w = 2000 kms" 1 and M w = 3.0 x 10~ 5 M yr _1 . 
This gives us /?st ~ 1.9 pc, well within the free-wind re- 
gion of the circumstellar bubble. This is the radius at which 
the production of the afterglow radiation is expected to end. 
In order to have a constant density medium during afterglow 
production, the wind termination shock has to lie much closer 
to the star, typically at ~ 0.1 pc (Chevalier et al. I2004i . If 
the entire afterglow shows a constant density medium, the jet 
has to pass the wind termination shock before it starts to de- 
celerate, i.e. before the end of the free expansion phase. The 
transition from free expansion to the relativistic Sedov-Taylor 
phase takes place at radius /?d (Eqn. \12i . With the values 
quoted above, and for 77 = 100, the deceleration radius equals 

= Rst/i 2 ~ 2 x 10~ 4 pc. In Section 3 we will consider var- 
ious possibilities for placing the wind termination shock this 
close to the star. 



2.3. Confining pressure 

The critical factor in bringing the wind termination shock 
closer to the star is the confining pressure of the surrounding 
medium. The pressure force constrains the expansion and size 
of the free-streaming wind region. The nature of this confin- 
ing pressure depends on the circumstances. For instance, in the 
scenario described in Sect. 12. H it is the thermal pressure in the 
bubble of shocked wind material, which in turn equals the ram 
pressure of the ISM associated with the overall expansion of 
the circumstellar bubble. For a rapidly moving star, it equals 
the ram pressure of the ISM associated with the stellar motion. 
One can generalize Eqn.|5]to: 



R _/MwVw\ I/2 



(17) 



with P c f the confining pressure. 

If we stipulate that the GRB jet has to pass the wind termi- 
nation shock before it enters the non-relativistic Sedov-Taylor 
expansion phase, this means that Rsj > Rv,- Combining 
Eas.ll5landll7lvields: 



Egrb > 



MwC 2 /M w V ff \" 2 



v w 



I 4/tPh-- ) 



Therefore, the confining pressure must satisfy: 



Pet > 



^ n ^GRB ^ w 



(18) 



(19) 



For the typical Wolf-Rayet wind values given in Sect. l2.2l and 
for £grb = 2 x 10 52 erg, this gives a confining pressure equal 
to: P cf =a 1.8 x 10~ 9 dyne cirT 2 . 

If, on the other hand, we demand that the jet has to pass 
through the wind termination shock before it has started to de- 
celerate, we get Ra > R w . This would have the advantage 
that the entire afterglow is generated in the constant density 
medium, but this situation is much more difficult to bring about. 
From the relationship between R$ T and R^ (Eqn. 1 1 51 we find 
that the confining pressure now must satisfy: 



4-n F 2 V 

GRB Vw 



(20) 



Since rj > 100, the confining pressure becomes extremely 
large: P c f =s 1 dyne cirT 2 for typical parameters. Such a large 
confining pressure is almost impossible to create. For instance, 
the thermal pressure of a static gas is equal to: 



P= 1.38xl0" 16 «r, 



(21) 



with n the particle density. So, for a particle density of 100000 
particles cm 2 , which is already extremely high, we would need 
a temperature of approximately 10 12 K. It is clearly not realistic 
to demand that the wind termination shock lies within the de- 
celeration radius. In a similar fashion, if the confining pressure 
is a ram pressure, 



P = p v 2 = 1.67 x 10" 



V 



1000 km/s 



(22) 
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GRB 

name 



redshift 



(x 10 52 erg) 



density profile 



GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 



970228 

970508 

980703 

990123 

990510 

991208 

991216 

00030 1C 

000418 

000926 

010222 

011121 

020405 

021004 

021211 



0.695 
0.835 
0.966 
1.60 
1.619 
0.706 
1.02 
2.03 
1.118 
2.066 
1.477 
0.36 
0.69 
2.3 
1.01 



2.68 <5) 
0.985 (5) 
8.09 <5) 
232 (5) 
21.9 (5) 
18.6 (5) 
7.83 (5) 
6.35< 5 » 
11.6 <5) 
45.6 (5) 
154< 6 > 

H <6) 

7.37 <7) 
2.2 (8) 
0.61< 6) 



W (D 

w (1.3.4) 
u <2) 
c (l,2,3,4) 
c ( 1,2.3,4) 
u (2) > w (3,4) 

W< 3 > 
u (3.4) 

W0. 4 ) 



v<3) ,,(4) 



Table 1. The redshift end isotropic energy for a number of 
GRBs as well as the density profile (w=wind, c=constant, 
u=undecided) of the circumstellar medium. If the density pro- 
file is marked as undecided, this means that both a wind and 
constant density profile are possible. Sources for the den- 
sity profiles are:(l) Chevalier & Li 2000, (2) Panaitescu & 
Kumar |2"0"0"T1 (3) Panaitescu & Kumar |2"0"0"21 and (4) Chevalier 
et al 120041 The redshifts were taken from the table at: 
http://www.mpe.mpg.de/~jcg/ maintained by J. Greiner. The 
values for the isotropic energy come from (5) Bloom et al. 
I200T1 (6) A mati l200ll (7) Price et al. l2003l and (8) Holland 
et al.ExE] 



it would require a large density in combination with a large 
velocity in order to achieve such a large value. Therefore, we 
will limit our discussion to the weaker condition dl9> . placing 
the termination shock radius well within the radius for non- 
relativistic Sedov-Taylor expansion. 

2.4. Observational indications 

The distance a GRB blast wave or jet will travel before it starts 
to decelerate depends on its initial energy £grb, and on the 
density of the surrounding matter. This means that a progenitor 
star that leaves a low-density wind, e.g. due to low metallicity, 
and which produces a high-energy GRB, might have a higher 
chance to show the GRB afterglow light curve for a constant 
density medium, since the jet can travel far before it starts to 
decelerate (Ean. ll2t . 

We have compiled a list of those GRBs for which in- 
dications about the circumstellar density profile of the CSM 
have been derived from the afterglow emission (Table Fig. 
13 shows the location of these bursts in the isotropic energy- 
redshift plane. We find the bursts with clear indications of a 
constant density medium preferentially at high isotropic energy 
and high redshift, while the bursts with a free wind type CSM 
favor low energies and redshifts. 

This figure indicates that afterglows seem to follow a trend 
consistent with the considerations of this section. Furthermore, 
there seem to be more afterglows that show a wind-like CSM 
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Fig. 2. Distribution of GRBs over the redshift-isotropic energy 
space (from TableQ. GRBs whose afterglow shows a constant 
density profile (circle) tend to have high energy and occur at 
larger redshifts than the afterglows which are formed in a free- 
streaming wind (square). The ones marked as triangles are am- 
biguous, which means that either different groups disagree on 
the nature of its afterglow, or no group has decided on its na- 
ture. If the symbol is large, several groups have reached the 
same conclusion. References for these afterglow models are: 
Chevalier & Li ( 



( 2000 1, Panaitescu & Kumar ( 200 1 ), Panaitescu 
& Kumar (120021 and Chevalier et al. (12004b . 



than afterglows with a constant-density environment. However, 
if afterglows with a constant density profile indeed occur on av- 
erage at higher redshifts, the relative numbers of the two groups 
may change with the advent of new data from the SWIFT satel- 
lite. Also, in various individual cases, which are not considered 
here, a constant density profile is preferred (Piro et al. 2005 ). 

If the shocked wind does indeed provide the constant den- 
sity medium observed in GRB afterglows, one would expect to 
see observational evidence that the GRB jet crosses the wind 
termination shock, as may have been observed in the afterglow 



of GRB 050904 (Gendre et al.|2005). However, the transition 
from the free-streaming wind to the constant density medium 
may not be clearly visible in the afterglow lightcurve (Nakar & 
Granot. l2006l . 



3. Ways to get the wind termination shock close to 
the star 



The results derived in Sect. 12.11 provide us with a variety of 
ways in which the wind termination shock can be brought 
closer to the star. We have tested several scenarios, either an- 
alytically, numerically or both. For the hydrodynamical simu- 
lations we have used the ZEUS 3D code of Stone & Norman 
JT992l >. 
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3.1. Low metallicity 

The winds of hot massive stars depends on their metallicity. For 
Gamma-Ray Burst progenitors, the mass loss rates might be 
lower than for Galactic Wolf-Rayet stars, since they may occur 
at low metallicity (Wijers 2001 Langer & Norman 2006 ). Such 
stars tend to have lower mass loss rates during their Wolf-Rayet 
phase (Crowther et al. l2002l Vink & de Koter l2005> . Since the 
ram pressure of the wind depends linearly on the mass loss rate, 
this helps to bring the wind termination shock closer to the star. 

From Eqn.[5]we know that a change in mass loss rate of two 
orders of magnitude brings the wind termination shock only 
one order of magnitude closer to the star. Following Vink & 
de Koter J2005> . Yoon & Langer d2005t and Yoon, Langer & 
Norman ( 2006 ) give a metallicity dependence of the mass loss 
rate for Wolf-Rayet stars as: 



(23) 



This implies that, in order to go from 10 pc to 0.1 pc, 
the mass loss rate would have to be four orders of magnitude 
lower. While it is not likely that many observed bursts have 
such low metallicities, values as low as Z o /100 have been sug- 
gested based on observations (Chen et al. 2005 1. 

In principle, one has to consider that the metallicity depen- 
dence of the mass loss rate applies to the main sequence phase 
of GRB progenitors as well. A star with a weak Wolf-Rayet 
wind had a weak main sequence wind before. This means that 
the thermal pressure in the hot bubble will be low (Eqn. [3}, 
which puts the wind termination shock further away from the 
star. This effect decreases the overall effect of the weak Wolf- 
Rayet winds, but it is not strong: the inner radius of the wind 
bubble depends only on main sequence mechanical wind lu- 
minosity to the power of 0.2 (Eqn. 5). The dependence of the 
radius of the wind termination shock on the metallicity of the 
progenitor star was investigated by Eldridge et al. (j2p05l. For 
stars at lower metallicity, the wind termination shock tends to 
lie closer to the star,but the dependence is weak. 

A benefit of low mass loss rates is the low density of the 
wind, which conforms to many of the afterglow lightcurve 
models (Chevalier et al. 20041. This allows the GRB jet to pen- 
etrate further into the CSM before it decelerates (Eqn.ll5>. This 
can explain why a signature of the wind termination shock does 
not appear in many GRB afterglow lightcurves: the jet is al- 
ready past this point when the afterglow starts. 

3.2. A high-density ISM 

If the density of the ISM is increased, the whole wind-blown 
bubble becomes smaller, including the radius of the wind ter- 
mination shock. A large range of densities in the ISM surround- 
ing GRB progenitors is to be expected. Stars are formed inside 
molecular clouds, which means that the star may find itself in a 
high density region. Unfortunately, the dependence of the wind 
termination shock radius on the ISM density po is not strong, 
R w oc p 3 ^ 10 (c.f. Eqn. [5}- This means that the density has to be 
very high in order to bring the wind termination shock as close 
to the star as we need, see for instance Eldridge et al. (2005 1. 




20 

R[pc] 



Fig. 3. The density of the gas around a 35 M star with the 
density of the ISM set to po = 10~ 225 g crrT 3 (upper segment) 



and po = 10 gem (lower segment). The high density of 
the ISM in the second segment limits the size of the bubble 
considerably, which increases the density in the bubble. Both 
segments show the same epoch, just before the star explodes as 
a supernova. 




Fig. 4. A blow-up of the central portion of Fig. [3] In the upper 
segment the wind termination shock can be seen at R ~ 14 pc. 
In the lower segment the termination shock is located at R ~ 1 
pc. 
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Table 2. Assumed parameters for the three evolutionary stages of our 35 M model (Z=0.20). 
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Phase 


End of phase [yr] 


M w [Moyr 1 ] 


V w [kms- 1 ] 


^photon [S ] 


Main Sequence 


4.745x10 s 


6.6xl0- 7 


2 671 


3.4xl0 47 


Red supergiant 


5.098xl0 6 


5.5xl0- 5 


15 


3.3xl0 41 


Wolf-Rayet 


5.295x10 s 


2.8xl0~ 5 


2140 


3.3xl0 47 



Figures[3]and|4]show the effect of a massive increase in the 
density of the ISM. Both segments of the figure show the den- 
sity around a 35 M star based on the star models calculated by 
Schaller et al. ( 1992 ). In one segment, the density has been set 
at 1CT 22 5 gem -3 ; in the other it is 1CT 195 gctrT 3 . Both sim- 
ulations were done using the method described in van Marie 
et al. ( I2004ll2005all2005b> . For the simulation with the low- 
density ISM the resolution is equal to 1000 radial grid points 
by 200 angular grid points. The high-density ISM simulation 
has a grid of 500 radial grid points by 200 angular grid points. 
The main-sequence and the start of the red supergiant phase are 
simulated in ID, and the results are mapped onto the 2D grid. 
The end of the red supergiant phase, as well as the Wolf-Rayet 
phase, are simulated in two dimensions. 

The wind parameters for the three phases of the stars evo- 
lution can be found in Table [2 As can be expected from Eqn. 
13 a change in the ISM density by three orders of magnitude 
changes the location of the wind termination shock by approxi- 
mately a factor of 10. To bring the wind termination shock even 
closer to the star would be difficult. To get a termination shock 
distance R w ~ 0.1 pc one needs an ISM with a density of about 
po ~ 10~ 16 gem -3 , which is rather unlikely. Again, as in Sect. 
13.11 this scenario can only account for a few GRB progenitors. 



3.3. A high-pressure ISM 

The analytical solution described in Sect. 12. H is only valid if the 
ISM is cold so that the circumstellar bubble expands superson- 
ically with respect to the surrounding ISM. For a sufficiently 
high ISM temperature, the thermal pressure of the ISM will act 
as an extra confining force on the expanding shell. An exam- 
ple of such a situation is a circumstellar H II region. If the star 
photo-ionizes the surrounding gas, the shell will expand more 
slowly. The expansion speed may even become subsonic. In 
that case the shell dissipates, and the circumstellar bubble will 
have a density discontinuity at the location where the shocked 
wind material meets the photo-ionized ISM. This effect has 
been described in van Marie et al. (120041 12005a1l2005b> . 

A hot ISM is the explanation offered by Chevalier et al. 
( 2004 1 for the presence of matter with constant density close to 
the progenitor star. They describe the structure of a circumstel- 
lar bubble that has expanded into a low-density, extremely hot 
medium. In such a situation, the main sequence wind can never 
create a supersonically expanding shell. The expansion of the 
bubble remains subsonic, which means that the thermal pres- 
sure in the hot wind bubble must approximately equal the ther- 
mal pressure in the ISM. Pressure balance between the wind 
ram pressure and the pressure in the ISM, p w V 2 ~ Po, gives 
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Fig. 5. The mass density (solid curve) and the internal energy 
density (dotted curve) of a circumstellar bubble that is expand- 
ing into a high-pressure ISM. Shown is the situation at a mo- 
ment just before the star explodes as a supernova. Moving out- 
ward from the star (from left to right) we see the free-streaming 
Wolf-Rayet wind, the wind termination shock (at /? w ~ 1 pc), 
the hot bubble filled with the shocked Wolf-Rayet wind mate- 
rial, the Wolf-Rayet wind driven shell (at R ~ 5 pc), the bubble 
of shocked main sequence wind material, and the ISM (R > 6.5 
pc). 



log(d) [g cm 3 ]= -26.0 -25.1 -24.1 -23.2 -22.2 -21.3 -20.3 



0.25 




0.2 0.3 

X[pc] 

Fig. 6. The density of the gas around two massive stars dur- 
ing the main sequence. The CSM bubbles of the two stars have 
merged, creating a superbubble. The image corresponds to the 
situation approximately 800 years after the start of the evolu- 
tion of the bubble. 
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Fig. 7. The internal energy density of the CSM surrounding the 
two stars at the same epoch as in Fig. |6] The internal energy 
density is highest between the stars, where the stellar winds 
from both stars contribute to the heating of the gas. 



the radius of the wind termination shock as: 

1/2 



AttPq 



(24) 



0.1 cm 3 , close 
20041. This situation 



We have performed a simulation of a similar situation, the re- 
sults of which are shown in Fig. |5] We used the same 35 M Q 
star model that was used for the high-density simulation that 
was described in Sect. 13.21 The temperature of the ISM was 
set to 10 8 K. The ID radial grid has 10000 grid points. In Fig. 
[5] the wind termination shock lies quite close to the star, once 
again at R w ~ 1 pc. The Wolf-Rayet wind-driven shell is still 
visible, even at the moment of supernova explosion. The reason 
is the high density and pressure of the main-sequence bubble. 
The pressure inside the shocked wind is high, since it has to 
be equal to the ISM pressure. The density in this region is high 
because the whole bubble is very small: pb ~ 10~ 25 gcirT 3 
which corresponds to a particle density «h 
to the value quoted by Chevalier et al 
results in a slow expansion speed of the Wolf-Rayet shell. In 
a 2D or 3D simulation, the Wolf-Rayet shell would have been 
fractured by its collision with the Red Supergiant shell (Garcia- 
Segura et al. !1996a|[T996bl van Marie et al. l2005all2003b1 left 
by the progenitor in an earlier evolutionary phase. 

However, a circumstellar environment such as described 
here is difficult to create. Even if we consider a large associ- 
ation of massive stars that creates a superbubble, the merged 
wind-blown bubbles of a number of stars, it is almost impossi- 
ble to increase the pressure to the level required to push the 
wind termination shock back to R w ~ 0.1 pc. Equation 1241 
shows the problem quite clearly. If the wind termination shock 
has to be ~ 10 times closer to the star, the pressure in the sur- 
rounding ISM as to be a factor of ~ 100 larger. This means that 
the total mechanical luminosity of the association that feeds the 
bubble has to be ~ 10 5 times larger than the typical mechan- 
ical luminosity of a single star, as the pressure in the bubble 
scales with the total mechanical luminosity as Pb K c.f. 
Eqn.|3 Effectively, this means that 100,000 massive stars have 



to release the kinetic energy of their wind into the bubble si- 
multaneously. 

The only place where the gas pressure of the interstellar 
gas can reach the required level is in the interior of starburst re- 
gions. There, stellar winds and supernova explosions add suffi- 
cient energy to the surrounding gas to reach these temperatures. 
It should be noted, though, that the wind termination shock in 
Fig.[5]is still more than one parsec away from the star. To bring 
it closer by another order of magnitude would mean that the 
pressure of the surrounding gas has to be increased by two or- 
ders of magnitude. 

3.4. Colliding winds 

In reality the problem shown in Sect. 13.31 may not be insur- 
mountable provided two stars are sufficiently close together, for 
instance in a binary. In that case, the winds of the two stars col- 
lide directly, rather than simply feeding energy into the same 
superbubble. In the combined wind bubble the temperatures 
will rise very rapidly in regions where the winds collide. We 
will consider the binary scenario, in the following paragraph. 
In this paragraph we will consider the simpler case of two sta- 
tionary stars, which are not gravitationally bound. 

Figs. [6] and show what happens when two stars inde- 
pendently feed mechanical energy into the same bubble. The 
simulation was done in 2D on a cylindrical grid of 1000x500 
grid points. The two stars both have the same mass loss rate 
(M = 4x 10~ 6 Mq yr 1 ) and wind velocity (V w = 2000 km s" 1 ). 
The two stellar winds have been simulated by filling two half 
spheres on the Z-axis with gas, moving away from central 
points (the stars) with a stellar wind density distribution (de- 
creasing with the radius squared). The simulation includes ra- 
diative cooling, but the effect of photo-ionization is not con- 
sidered. The winds never collide directly. Instead, each star 
creates its own bubble. When their respective main sequence 
shells hit each other, the bubbles merge. The region directly 
between the two stars is unrepresentative of the rest of the re- 
gion, since the pressure is higher. However, the difference is not 
enough to bring the termination shock close enough to the pro- 
genitor star, unless the GRB is aimed directly at the second star. 
Still, if this were the case, the afterglow would not show a con- 
stant density medium, since the GRB jet would pass through 
the shocked wind material and run into the free-streaming wind 
of the second star. Should the winds of several stars collide di- 
rectly, then the pressure in the bubble would locally be very 
high, which would indeed constrain the radius of the wind ter- 
mination shock. The problem remains that the GRB jet would 
pass quickly through the layers of shocked wind material and 
the afterglow lightcurve would most likely show a variety of 
density profiles. 



3.5. A binary companion 

The presence of a binary companion can bring the wind ter- 
mination shock closer to the star For instance, in a Wolf- 
Rayet/O star binary system, the Wolf-Rayet star, which has the 
smaller mass, will orbit the O-star. The Wolf-Rayet wind will 
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Table 3. Wind parameters for the binary stars simulation Table 4. Assumed parameters for the stellar wind in our bow- 



shown in Fig.[S] 



shock model ('Figs.l9lto ll li . 





M w [Moyr 1 ] 


V w [kms- 1 ] 


phase 


M w [Moyr 1 ] 


K w [kms- 1 ] 


Star 1 


l.OxKT 5 


1000 


Red Supergiant 


l.OxlO- 4 


10 


Star 2 


l.OxlCT 4 


10 


Wolf-Rayet 


l.OxlO^ 5 


1000 












log(d) [g cm 3 ] s -28.0 -26.4 -24.B -23.2 -21.7 -20.1 -18.5 



X [AU] 



Fig. 8. The density of the gas around two stars with different 
wind parameters. The left star (at X = 7 AU) has a high 
mass loss rate and a low wind velocity, while the right star (at 
X — 9 AU) has a low mass loss rate but a high wind velocity, so 
that its ram pressure is the largest. The stronger wind (from the 
star on the right) dominates the CSM, while the weaker wind 
(on the left) is folded back, away from the companion star. 



be distributed over a larger volume, which decreases the den- 
sity of the wind and its ram pressure. However, this effect is 
probably not very significant. If we are dealing with a close bi- 
nary, the radius of the orbit is smaller than R w and the increase 
in volume is small. If we are dealing with a wide binary, the 
orbital motion will be so slow that the effect on the wind is no 
longer significant: to lowest order the Wolf-Rayet star can be 
treated as motionless. 

The wind emanating from the companion star is unlikely 
to be of any substantial influence: the ram pressure of a Wolf- 
Rayet star wind is much larger than that of the wind from nearly 
any other star. The wind from the companion will be pushed 
back in a bowshock, and its material will be constrained to a 
narrow layer in the orbital plane. This is demonstrated in Fig. 
[S] This is a simplified model of two stars in a binary system. 
The star on the left has a high mass loss rate, and a low wind 
velocity, like a red supergiant or luminous blue variable; the 
star on the right has a lower mass loss rate, but a high wind ve- 
locity, like a Wolf-Rayet star (Table QJ. The simulation is done 
in two dimensions, assuming a cylindrically symmetric situa- 
tion. At the beginning of the simulation the grid is filled with 
ambient material at a constant density. The simulation includes 
radiative cooling, but the effect of photo-ionization is not con- 
sidered. This simulation does not take the orbital motion of the 
stars into account. It shows quite clearly how the weaker wind 
is being folded back by the strong Wolf-Rayet wind. An ana- 
lytical solution pertaining to this situation was given by Canto 




Fig. 9. The density of the gas around a moving red supergiant 
star. The velocity field in the ambient medium is shown as 
stream traces. The bowshock shell shows Kevin-Helmholtz in- 
stabilities. 



et al. ( 1996 1, while numerical simulations have been presented 
by Stevens et al. (Q3532). 

The jet associated with the GRB is expected to move along 
the rotation axis of the progenitor star, which should be per- 
pendicular to the orbital plane. Therefore, the companion star 
and its wind would have little or no influence on the GRB jet 
and the afterglow light curve. Of course, it is possible that a jet 
propagates through the wind of another star, which happens to 
be nearby. In that case, the density profile encountered by the 
GRB would depend entirely on the angle at which the jet hits 
the free-streaming wind. If it would meet the companion wind 
head-on, it would actually encounter a density that increases 
as the jet penetrates deeper into the wind. 

3.6. Stellar motion 

If the progenitor star is moving rapidly through the interstellar 
medium, the situation becomes completely different. The equa- 
tions in Sect. 12. II apply when the star is motionless relative to 
the surrounding gas. However, if the star is moving at super- 
sonic velocity with respect to the ISM, the collision between 
stellar wind and ISM will create a bowshock. This bowshock 
forms at the distance from the star R^, where the ram pressure 
of the wind equals the ram pressure of the ISM in the frame 
where the star is at rest. The stand-off distance of the bowshock 
equals 



MwVw 
4np v\ ' 



(25) 
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Fig. 10. Similar to Fig.|9j but after the star has become a Wolf- 
Rayet star. The newly formed bowshock-shell system, which is 
sweeping up its predecessor, is highly unstable. 
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Fig. 11. The internal energy density at the same moment in time 
as in Fig. ^| The wind termination shock is clearly visible, 
0.5... 1.25 pc away from the star, depending on the angle. 



where v* is the velocity of the star (Chen et al. 1995 and Wilkin 
1996 1. For a Wolf-Rayet wind, with a velocity of 2 000 km s _1 
and mass loss rate 3x 10~ 5 M Q yr _1 , from a star that moves with 
a velocity of 80 km s _1 through an ambient medium with a den- 
sity of 10 _22i gem -3 , this gives us a radius of the bowshock 
directly in front of the star of R\, s =s 1 .25 pc. Generally, the 
radius of the bowshock depends on the angle between the po- 
sition vector and the direction of motion. Similarly, the radius 
of the reverse shock is angle dependent, since the hot bubble 
of shocked wind material is not isobaric. The thermal pressure 
will be highest directly in front of the star, where the wind col- 
lides head-on with the ambient medium. 

Such a rapidly moving star can be the survivor of a close 
binary, where the companion has exploded as a supernova. In 
this example, the force of the supernova explosion broke up 
the binary, and the surviving star has moved away at approxi- 
mately its original orbital velocity. Such runaway stars are not 
uncommon, but they do not always form a bowshock. If the 
star moves through a hot region of space (such as an H II re- 
gion, or the circumstellar bubble of another stars), its motion is 



subsonic and no bowshock is formed (Huthoff & Kaper l2002l . 
During the early stages of its runaway phase, the star cannot 
form a bowshock, since it is moving through its own main- 
sequence wind bubble, where the temperature will be of the 
order of 10 7 K, and the speed of sound is of the order of 300 
km/s. How much time the star needs to escape its own bubble is 
extremely difficult to predict, as both the velocity of the star and 
the radius of the bubble can vary by at least an order of mag- 
nitude. Typically, this will be of the order of 10 5 ...10 6 years. 
Several Wolf-Rayet stars with bowshocks have been observed, 
such as Ml -67 around WR 124 (van der Sluys & Lamers l2003l 
and possibly NGC 2359 around HD 56925 (Chen et al. lT995V . 

To better understand the morphology of the circumstellar 
medium around a fast moving star, we have performed a two- 
dimensional hydrodynamical simulation. Since the stellar wind 
is assumed to be equally strong in all directions, it is possible 
to use a cylindrical grid, using the <£ axis as the symmetry axis. 
The Z-axis coincides with the direction of motion. This is sim- 
ilar to the approach used by Brighenti & D'Ercole ( 1997 1 and 
Comeron & Kaper (1998 1. Initially, the grid is filled with an 
ambient medium (po = 10~ 225 g cm -3 ), which travels along 
the Z-axis at a constant velocity (v* = 80kms -1 ). Matter 
flows in at one end of the grid, and out at the other. The wind 
from the star is simulated by filling half a sphere on the Z-axis 
with material, which moves outward from a central point, and 
which has the density distribution of a stellar wind (as we did to 
model the binary system in Sect. l3.5> . The grid resolution has to 
be rather high, in order to resolve the Kevin-Helmholtz instabil- 
ities in the shell (Brighenti & D'Ercole ll997> . The stellar wind 
in our models follows two separate evolutionary phases: first 
the red supergiant phase, followed by a Wolf-Rayet phase. The 
wind parameters for these two phases can be found in Table |4] 
Photo-ionization is not included in this simulation. 

In Figs. [9] through Fig. ^2 we show the results. Figure [9] 
displays the density of the circumstellar medium during the 
red supergiant phase. The bowshock shell is clearly visible: it 
shows a series of so-called 'Cat's eyes', which are the result of 
a fully-developed Kevin-Helmholtz instability. These instabili- 
ties only occur if the ambient medium is cold. If the star moves 
through an HII region, the instabilities disappear. In Figs. 1101 
and[^th e star and wind are in the Wolf-Rayet phase. The fast 
Wolf-Rayet wind has created a new shell, which sweeps up the 
old shell. It is extremely unstable. In Fig.[^ which shows the 
internal energy density of the CSM, we can see the location of 
the wind termination shock. The pressure in the CSM is highest 
in front of the moving star, where the collision between wind 
and ambient medium is the most violent. The radius of the wind 
termination shock varies between 0.5 and 1.25 pc. 

Typically, one would expect the star to move in what was 
originally the orbital plane. This means that the GRB would be 
directed about 90" away from the direction of motion, since the 
jets should move along the stellar rotation axis. The radius of 
the wind termination shock at that position angle is about 1 pc. 

4. Combining different scenarios 

It is possible that several of the mechanisms considered in Sect. 
|3]act together. However, not all combinations are possible. 
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Fig. 12. The analytical solution for the radius of the reverse 
shock around a Wolf-Rayet star, depending on the parameters 
of the Wolf-Rayet wind and the density of the ISM. The radius 
of the reverse shock was calculated with Eqn. [5] The Wolf- 
Rayet lifetime has been taken as constant and the wind parame- 
ters for the main sequence and red supergiant phase come from 
Table |3 The parameter space in this figure covers Wolf-Rayet 
wind velocities between 500 and 5000 kms -1 , mass loss rates 
between 10~ 7 and 1.0~ 4 M yr~' and ISM densities ranging 
from 5 x 10~ 4 to5 x 10 4 particles per cirT 3 . The iso-surfaces 
in the plot show the planes with a constant radius of the wind 
termination shock. 



4.1. Low metallicity and high ISM density 

A low ram pressure of the wind can be combined with any other 
mechanism, since it depends only on the internal structure of 
the star, while the other mechanisms rely on special external 
conditions such as the density or pressure in the ISM. The ef- 
fect of combining different wind parameters and variations in 
the density of the ISM is shown in Fig. El Here we have taken 
a massive star and varied the wind velocity and mass loss rate 
during the Wolf-Rayet phase. As a basis we use the 35 M Q 
star model from Table |2] The evolutionary time and the pa- 
rameters for main-sequence and red supergiant wind are kept 
constant, but the parameters for the Wolf-Rayet wind are var- 
ied. The wind velocity of Wolf-Rayet stars varies between 500 
and 5000 km s _I , while the mass loss rate varies between 1 .0~ 7 
and 1.0~ 4 M yr~' (c.f. Hamann et al. 1 19951 and Koesterke & 
Hamann 1995 1. We have plotted the mass loss rate, the wind 
velocity and the density of the ambient medium in Fig. [2] The 
density of the surrounding ISM was varied between 5 x 10~ 4 
and 5 x 10 4 particles per cm -3 . 

The results show a considerable variation in the location of 
the wind termination shock. As discussed above, the constant 
density inferred from some GRB afterglows is thought to be a 
signature of the shocked wind material. The quantitative value 
of this density can be calculated analytically, by dividing the 
total mass loss of the star over the spherical volume of the bub- 
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Fig. 13. The analytical solution for the radius of the bowshock 
in front of a Wolf-Rayet star, depending on the stellar wind pa- 
rameters, the velocity of the star and the density of the ISM. 
On the horizontal axis we give the parameters of the stellar 
wind, on the vertical axis, the ram pressure of the ISM. The 
radius of the bowshock was calculated with Eqn. [25] The wind 
velocity, mass loss rate and ISM density cover the same inter- 
vals as in Fig. [21 The stellar velocity varies between 10 and 
100 kms _I . The circle shows the location of the numerically 
calculated model from Figs.l9ltoll II 



ble. The outer edge of the bubble has a radius (Weaver et al. 
IT9T71 



R = 



/ 250L W \ 



1/5 



,3/5 



(26) 



Typically, the resulting particle density is less than 10 cm 3 
(see also Sect. 13.21 and the Figs. [5] and|4j, which is acceptable 



according to the observations (Stratta et al. 2004). 

4.2. Stellar motion in high density, and low metallicity 

In Fig.[D]we show how the radius of the bowshock varies with 
the stellar velocity v+ and the density of the interstellar medium 
po in the combination pqv\ (the ISM ram pressure), and the 
wind parameters of the star in the combination MV W , the abso- 
lute momentum discharge of the wind. The radius of the bow- 
shock follows from Eqn. [25] This parameter study is simpler 
than the case of the stationary star, since the wind velocity and 
mass loss rate do not have to be plotted independently. The 
variations in wind parameters and ISM density are the same as 
those used to produce Fig. [21 while the stellar velocity varies 
between 10 kms~' and 100 kms -1 . Fig.ll3lshows that the ra- 
dius of the bowshock changes considerably, depending on these 
parameters. The bowshock can be located closer to the star than 
the wind termination shock for a stationary star (see Fig.ll2i. 
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In this figure we also show the location of the model that 
we simulated in the manner described in Sect. 13.61 The numer- 
ical result from the simulation is quite close to the analytical 
approximation. The dependence of the bowshock radius on the 
ISM density po is much stronger than the density dependence 
of the wind termination shock radius for a stationary star, c.f. 
Eqns.|5]and|25j which makes this combination particularly ef- 
fective. 

The main problem with the results of Fig.^]is the fact that, 
for extremely small bow-shock radii, the constant density is 
provided by the ISM, rather than by the shocked wind material. 
In those cases the density in the ISM is too high compared to 
the observations. These place a limit on the particle density 
of about 10 ctrr 3 (Stratta et al. 2004). This problem does not 
occur for GRB afterglows produced in a low-density shocked 
wind. 

4.3. Other combinations 

Combining a high density with a high pressure in the ISM is 
nearly impossible, since a medium with a high thermal pressure 
will expand quickly (on its own hydrodynamical timescale), 
and the density will decrease rapidly. Actually, the high den- 
sity would have little or no effect. If the pressure is as high as 
envisioned in Sect. 13.31 then the thermal pressure of the gas is 
much stronger than the ram pressure of a moving shell. This 
means that the thermal pressure is the only relevant confining 
force that acts on the shell, and the actual density of the ISM is 
no longer relevant. 

Binary companions can occur in any environment, irre- 
spective of the density or pressure of the surrounding ISM. 
However, since the influence of a companion star on the loca- 
tion of the wind termination shock is limited at best, this makes 
little difference. 

Combining stellar motion with a high-pressure environ- 
ment does not lead to a viable scenario, as the velocity of the 
star would automatically become subsonic, and the motion of 
the star would no longer be a relevant factor. In any case, the 
runaway star would have to live long enough for the star to 
leave its own circumstellar bubble, since the temperature of the 
shocked wind material would again make the stellar motion 
subsonic. 

5. Conclusions 

Of the scenarios discussed, the most likely are: a Wolf-Rayet 
wind with low ram pressure, a high density ISM, and a mov- 
ing progenitor star. It is the last possibility that can bring the 
wind termination shock extremely close to the star, as shown 
in Fig.[0] Observations show that the wind velocity and mass 
loss rate vary considerably from one Wolf-Rayet star to an- 
other. Since massive stars are thought to form within molecular 
clouds, a high density ISM is also possible. Fast moving Wolf- 
Rayet stars have been observed inside our galaxy. The combi- 
nation of a weak Wolf-Rayet wind and a high ISM density can 
by itself be enough to obtain the desired result, but the choice 
of both parameters will have to be rather extreme. 



A high pressure ISM can produce the desired result as well, 
and probably accounts for some of the observed GRB after- 
glows. However, regions in the ISM with a sufficiently high 
temperature are rare. 

While many massive stars have binary companions, this by 
itself will hardly affect the location of the wind termination 
shock, since the ram pressure associated with the wind from 
a Wolf-Rayet is much stronger than that of any possible com- 
panion. The Wolf-Rayet wind will therefore dominate the mass 
flow from the binary system. Only if the binary companion 
were also a Wolf-Rayet star might its wind be competitive. 

Only one such system has ever been found: WR 20a (Rauw 
et al. 20041, which contains two WN type stars. However, since 
these two stars still have most of their hydrogen envelope, they 
are not yet GRB progenitors. Unlike most WN stars, these are 
main sequence stars with a very high mass loss rate. The orbital 
motion of a star in a binary system changes the properties of the 
stellar wind, but its effect is not enough to make a significant 
difference to the location of the wind termination shock. 

There is one possibility that has not been discussed so far: 
the constant density medium may not be in the hot bubble 
that drives the main sequence shell, but in the small bubble 
of shocked Wolf-Rayet wind material, which drives the Wolf- 
Rayet shell into the earlier Red Supergiant or Luminous Blue 
variable wind. For this to be true, the GRB would have to oc- 
cur very early during the Wolf-Rayet phase. (Typically, a Wolf- 
Rayet wind driven shell moves at velocities of the order of 100 
kms" 1 , so it reaches 1 pc within 10000 years.) Apart from the 
extremely narrow time frame, there are several other problems 
with this scenario: During the early stages of its development, 
the layer of shocked Wolf-Rayet wind material is extremely 
thin (< 1 pc), so any afterglow generated in this medium would 
almost certainly show signs of either the transition between un- 
shocked and shocked Wolf-Rayet wind, the Wolf-Rayet wind 
driven shell, or both. Second, the layer of shocked Wolf-Rayet 
wind, due to its thinness, rarely has a truly constant density, 
since both Rayleigh- Taylor instabilities of the Wolf-Rayet shell 
as well as internal turbulence in the shocked wind contribute to 
local density changes. 

The main problem of the scenarios discussed to bring the 
wind termination shock close to the star is that the parame- 
ters of stellar mass loss or the conditions in the ISM have to 
be pushed to extremes. Therefore, none of the mechanisms de- 
scribed here seems, by itself, to account for the high percentage 
(~ 25 %) of the observed characteristics of GRB afterglows. 
However, the large number of possibilities and, in particular, 
possible combinations of the scenarios presented here, may ex- 
plain why a significant fraction of GRB afterglows seem to oc- 
cur in a constant-density medium. 

Generally speaking, outside influences like the motion of 
the star and the condition of the ISM can bring the termination 
shock to about 1 pc. In order to bring it any closer, the Wolf- 
Rayet wind has to be weak, which means that the progenitor 
star most likely had low metallicity. From our models of the cir- 
cumstellar medium and the GRB afterglow models discussed in 
Sect. I2.4l we can conclude that the GRB afterglows generated 
in a constant density medium come from low metallicity stars 
with relatively strong GRBs (although it is an open question, 
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as to whether the GRB itself was strong or the GRB jet had a 
narrow opening angle). 
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